Figure 10 Distributions of rest equivalent widths in the LCDM model with J_ 2 i = 3.0 
at redshifts 2.0, 2.4, 2.8, 3.2, 3.6 and 4.0. The dashed line is given by Eq.(15) with 
W* = 0.27 A. 

Figure 11 Two-point correlation functions of Lya forest lines with W > 0.16A in veloc- 
ity space for model of LCDM with J_ 2 i = 3.0. The error bars represent 1 a. 

Figure 12 The Gunn-Peterson optical depth in the model of CHDM with J_ 2 i = 0.3, 

Figure 13 Number density dN/dz of Lya lines at W t h = 0.32 and 0.16 A as a function 
of redshift in CHDM model. The threshold S t h = 4.5. For dotted and dot-dashed 
lines, the UV background is J_ 2 i = 0.3, and the normalization of the amplitude 
of initial perturbation is taken to be Q rms _ps = 18 and 13fiK, respectively. For 
dashed and solid lines, the UV background is J_ 2i = 1.0, and the normalization 
of the amplitude of initial perturbation is taken to be Q rms _ps = 18 and 13fiK, 
respectively. The open circles are data from Bechtold (1994). 

Figure 14 Number density dN/dz and its ^-dependence in the models of LCDM with 
S t h = 4.5. The UV background is taken to be a) constant J_ 2i = 0.3 (long dashed 
line), and b) the model of QSO "turn on" at z = 5 (dot-dashed line). 
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Figure Captions 

Figure 1 Examples of simulated spectra of the Lya forest: (a) SCDM, J_ 2 i = 0.3, where 
J_ 21 is related to J c by J c = J_ 21 x 10" 21 cm" 2 s" 1 Hz" 1 sr" 1 ; (b) LCDM, J_ 21 = 1.0; 
and (c) CHDM, J_ 2i = 0.3. In all three spectra, the threshold parameter 8 t h = 1.69. 

Figure 2 Number density dN/dz of Lya lines at Wth = 0.32A as a function of redshift. 
The threshold 8 t h = 1-5. The UV background is J_ 2i = 0.3 (dotted line); 1.0 (dashed 
line); and 3.0 (long dashed line). The open circles are data from Bechtold (1994). 
Models (a), (b), and (c) are for SCDM, LCDM, and CHDM, respectively. 

Figure 3 Number density dN/dz of Lya lines at Wth = 0.16A as a function of redshift. 
The threshold 8 t h = 1-5. The UV background is J_ 2i = 0.3 (dotted line); 1.0 (dashed 
line); and 3.0 (long dashed line). The open circles are data from Bechtold (1994). 
(a), (b), and (c) are for models of SCDM, LCDM, and CHDM, respectively. 

Figure 4 The same as Fig. 2, but 8 t h = 4.5. 

Figure 5 The same as Fig. 3, but 8 t h = 4.5. 

Figure 6 The same as Fig. 2, but 8 t h = 7.5. 

Figure 7 The same as Fig. 3, but 8 t h = 7.5. 

Figure 8 Distributions of rest equivalent widths in LCDM with J_ 2i = 3.0. Simulated 
results: low redshifts 1.8 < z < 2.5 (dotted line); high redshifts 2.5 < z < 3.2 
(solid line). Observational data: z < 2.5 (open circle); z > 2.5 (solid circle). All 
distributions are normalized to the observed numbers at 0.32A. 

Figure 9 The same as Fig. 8, but for the model of CHDM. 
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does not significantly correlate to column density Njji and to redshift. However, some 
observations showed the redshift-dependence of b (e.g. Williger et al. 1994; Pettini et al. 
1990). 

Therefore, one cannot undoubtedly discriminate between SCDM and LCDM by the 
current simulation and real data. Nevertheless, a systematic difference between SCDM 
and LCDM can be seen significantly: good fitting of SCDM is generally given by higher 
threshold S t h than in LCDM. In other word, the SCDM model needs an undercutting of 
collapsed regions than the LCDM, i.e. it appears to favors a significant part of the Lya 
clouds to be related to halos of collapsed objects at high redshifts. The low redshift Lya 
absorption lines have already been found to be partially produced by galactic halos (e.g. 
Mo & Morris 1994). This problem now is also important in high redshifts. Because the 
number of the collapsed regions identified by, say, PS formalism should be equal to or 
larger than the number of observed collapsed objects at high redshift, such as CIV or 
Mgll metal systems in QSOs absorption spectrum, the occurrence of halo related Lya 
forest absorptions might be estimated by the number density of these metallicity systems 
(Bi & Fang 1995). 

We thank Dr. Y.P. Jing for useful discussions. JG thanks also Dr. J. Bechtold for 
her helps. HGB thanks World Lab for a fellowship. This work was partially support by 
the NSF INT-9301805. 
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This result is in a good agreement with other studies on larger scale structure forma- 
tion in the CHDM model. The "standard" CHDM model of Vt c /Vt h /Vt h = 0.6/0.3/0.1 was 
already found to be lack of perturbation powers to form high redshift collapsed objects. 
The N-body simulation of CHDM at moderate redshift showed that a rapid formation of 
collapsed objects appeared only in the near past z ~ 0.5 (Jing & Fang 1994). CHDM is 
also difficult in explaining the formation of damped Lya systems (Klypin et al. 1994), and 
high-redshift halos (Ma & Bertschinger, 1994). Therefore, our results strengthened the 
conclusion that the "standard" CHDM model, i.e. 60% cold and 30% hot dark matters 
and 10% baryons, would be in problem to pass tests of high redshift objects. 

When J 2 i < 0.3, all results of dN/dz curves in LCDM model are larger than real 
data. Considering the linear approximation probably underestimate the number of clouds 
at high redshift, the simulated dN/dz of LCDM would be a lower limit to the real number. 
When non-linear evolution is considered, the simulated result will be more larger than 
real data. Therefore, J 2 \ < 0.3 should be a constraint to the model of LCDM. 

The deficit of the W > 0.32A lines at z < 2.5 in J 21 = 3.0 LCDM (Figs. 2b, 4b and 
6b), and the deficit of W > 0.5 A objects in z < 2.5 width distribution of Fig. 8 can be 
addressed on the uncertainties of parameters. If we replace the constant UV background 
by an evolutionary model, say the "turned on" model of QSOs at redshift z = 5 (Bajtlik, 
Duncan & Ostriker 1988; Bechtold 1994), the deficit disappeared (see, dot-dashed lines 
in Fig. 14). Therefore, within a reasonable range of parameters, LCDM can reproduce all 
observational properties of the Lya forest. 

We should also pointed out the uncertainties in the current observations. First, the 
measurements of dN/dz, especially at z ~ 4, are still dispersive. Although the evolution of 
dN/dz given by different groups (Lu et al. 1991; Press & Rybicki 1993; and Bechtold 1994) 
showed a common trend, their values of dN/dz are differences by about 30%. Second, we 
have taken the velocity broadening parameter b to be constant. This assumption is based 
on the observation that the distribution of 6, well concentrated around 25~30 km s _1 , 
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from QSO 1202-0725 (z em = 4.7). Therefore, CHDM cannot pass the Gunn-Peterson test. 
This shows once again that the late structure formation in the CHDM model cannot be 
consistent with observations at high redshift. 

4. Conclusions and Discussions 

We simulated Lya forests in QSO's spectra in linear regime, in which the Lya clouds 
are modeled as the absorption of HI in less developed regions of the density field. Collapsed 
regions in the density field are identified, and then removed by a Press-Schecter-like 
criterion, i.e. a threshold of density contrast. The simulated samples have been compared 
with the observational data in four aspects: 1) the number density of Lya lines and its 
dependencies on redshift and equivalent width; 2) the distribution of equivalent widths 
and its redshift dependence; 3) clustering; and 4) the Gunn-Peterson effect. 

In spite of uncertainties in parameters of the UV background J 2 \ and the contrast 
threshold <5, our simulation has shown that the Lya forest can set effective constraints to 
models of large scale structure formation. The current observations have already given 
almost no room to the CHDM model even when the large uncertainties in the parameters 
are considered. In fact, one can see from Fig.l that the simulated spectrum in the CHDM 
model fails to show the forest-like features. The square of variance of the baryonic matter 
in CHDM model is < 0.5 in redshift range 2 < z < 4. Non-Linear evolution is negligible. 

This result still hold for a larger normalization of the amplitude of the density per- 
turbations. The two years of COBE-DMR observation (Bennett 1994) showed that the 
amplitude of the quadrupole anisotropy of cosmic background radiation might be as large 
as Q rms _ps = 18 ± 1.5 fiK. Using this amplitude we repeated the simulation. The results 
for the CHDM model are plotted in Fig. 13, which contains two sets of the curves dN/dz 
vs. (1 + z), corresponding to the normalization of Q rms _ps = 18 and 13fiK, respectively. 
One can see from Fig. 13 that the larger normalization do rise the number of Lya clouds. 
However, there are still much less than observations. 
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weakly each other, even the density field is produced by a power law spectrum of the 
perturbations. Since collapsed regions have been removed, one can expect that the line- 
line correlations would be more weak. 

The two point correlation functions of W > 0.16A lines for various redshifts in the 
LCDM model are plotted in Fig. 1 1 . In the entire redshift range from 2.0 to 4.0, no signifi- 
cant structures are detected. Therefore, the simulated spectra are, indeed, consistent with 
what observed, i.e. no cloud-cloud correlations exist on scales from 100 km s -1 to 2000 
km s -1 . It should be pointed out that the comparison cannot be taken seriously on scale 
less than 100 km s -1 or 1 h _1 Mpc, because of the line-blending caused by the velocity 
dispersion and the point-spread-functions of FWHM~ 30 km/s. As a consequence, any 
clustering on that scales will essentially be removed. An anti-correlation is then to appear 
in Fig. 1 1 . It should also be pointed out that the conclusion of no line-line correlation make 
sense only statistically. It does not exclude rare events, such as possible voids (Dobrzycki 
& Bechtold 1991). 

3.4 The Gunn-Peterson test 

The simulated density contrast 1 + (x, z) contains hydrogen gas in both clumpy and 
diffuse forms. Therefore, the Gunn-Peterson optical depth, which came from the diffuse 
component, can directly be estimated from the minimums in t(Z). No recognizable 
Gunn-Peterson depth can be found in the models of SCDM and LCDM. Assuming the 
signal-to-noise ratio S/N ~ 29 for the continuum in the simulation samples, the upper 
limit of the depth is ~ 0.01 from z = 1.6 to 4 in SCDM and LCDM. 

However, for the model of CHDM with J-21 = 0.3, which is relatively the best 
one among all fluxes being considered, the Gunn-Peterson optical depth is found to be 
tqp ~ 0.04 at z = 3.4; 0.07 at z = 3.6; and 0.10 at z = 4.0 (Fig. 12). From a measurement 
of t gp by spectrum of QSO 0000-263 (z em = 4.11), we have t gp < 0.04 at z w 3.8 (Webb 
et al 1992). A new limit, tqp < 0.02 at z ~ 4.3, is obtained by Giallongo et al. (1994) 
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al. 1986; Rauch et al. 1992; Weymann 1992). W* may also depend on redshift in the 
sense that the distribution steepens at lower redshift (Bechtold 1994). 

Fig. 8 shows the probability distributions of the rest equivalent widths in the J_ 2 i = 3 
and 6 t h = 4.5 LCDM model. The solid curve is the mean from the samples at z = 
2.6,2.8,3.0 and 3.2; the dotted curve from the samples at z = 1.8,2.0,2.2 and 2.4. The 
theoretical distribution n{W) has been normalized to the observation (Bechtold 1994) at 
W = 0.32A; this is equal to normalize the total observed W > 0.16A lines at z = 2.44. 
The agreement between the model and the observation is rather good. This result cannot 
simply be explained as a less dependence of the width distribution on models. Actually, 
the width distribution is sensitive to models. For instance, we calculated distribution of 
the equivalent widths in the model CHDM with J -21 = 0.3 and S t h (Fig. 9). Both the 
z > 2.5 and z < 2.5 distributions are much steeper than the observational data. There 
is almost no W > 0.6A lines produced at z < 2.5. Therefore, the CHDM model is also 
uncomfortable in the test of the width distribution. 

The W distributions at individual redshifts 2.0, 2.4, 2.8, 3.2, 3.6 and 4.0 in J 2 \ = 3 
LCDM model are plotted in Fig. 10. The dotted line in Fig. 10 is given by eq.(15) with 
W* = 0.27A. A rapid evolution of the width distribution is found, especially at era after 
z = 2.4. 

3.3 Clustering properties 

Lya clouds show an absence of two-point correlation in their spatial distribution. 
Some early samples were found to be excessive correlation at velocity differences less 
than ~ 300 km s -1 (Webb 1987). Later samples usually detected no clustering at all. 
This fact has been used as a negative grade of the model of dark matter confinement, 
because it was argued that if Lya clouds formed in the density field of dark matter, 
the spatial correlations of the clouds should be about the same as galaxies (e.g. Bajtlik 
1992). However, it has already been shown in B93 that the simulated Lya lines correlate 
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Therefore, the CHDM model is difficult to pass the test concerning subgalactic structures 
at high redshift. Actually, in the case of the CHDM model, the square of variance of the 
baryonic matter at the considered high redshift is only about 0.5, i.e. it is truly in the 
stage of linear evolution. Therefore, the failure of CHDM would not easily be overcame 
by non-linearity. 

The results of SCDM sensitively depend on S t h- Figs. 2a - 7a showed that the better 
fittings are given by S t h > 4, i.e. less removal of the collapsed objects. It implies that in 
the SCDM model there are too much powers of the density perturbation on scales of few 
to tens Mpc, and baryons were already dragged into collapsed regions at very early time. 
Therefore, the baryons remaining in pre-collapsed regions are not enough to contribute 
to Lya clouds. Therefore, the SCDM model requires a significant fraction of Lya forest 
lines produced by gas component located in the regions of collapsed objects. 

As has been known, comparing with the SCDM, the LCDM model has lower per- 
turbation powers on scales of few kpc to about 10 Mpc, an overproduction of collapsed 
objects could be avoided. It will give enough pre-collapsed baryons to form the Lya forest. 
On the other hand, unlike CHDM, LCDM predicts a significant development of clustering 
at z > 1, Lya forest systems and other absorption systems could be formed early enough. 
As a result, dN/dz in the LCDM model shows good fitting for all reasonable ranges of 
parameters (Figs. 2b - 7b). Both the W t h = 0.32A and 0.16A curves of dN/dz with 
J-2i = 3 in LCDM are consistent with the observations within 1 a. 

3.2 Distribution of equivalent widths 

Within a given redshift range, the distribution of the rest frame equivalent width W 
is given by 

n(W) cx e- w ' w * , (15) 
where W* is found to be from about 0.23 A to 0.35 A (Sargent et al. 1980; Murdoch et 
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as a function of z in the following discussion. In this paper, we are limited to study high 
redshift Lya forests (z > 1.6). Low redshift Lya absorption forests showed, sometimes, 
different features from that at high redshifts and might be questionable to be described 
by a linear theory. 

The simulated and observed ^-dependence of dN/dz are plotted in Figs. 2, 4 and 6 
for W > 0.32A, and Figs. 3, 5 and 7 for W > 0.16A. The observational data are taken 
from Bechtold (1994). The UV background is assumed to be constant with flux J-21 = 
0.3 (dotted), 1.0 (dashed) and 3.0 (long dashed). The threshold contrast S t h is taken to 
be 1.5 (Figs. 2, 3), 4.5 (Figs. 4, 5) and 7.5 (Figs. 6, 7). 

Obviously, comparing with the PS formalism, S t h > 3 would be too large for identify- 
ing collapsed regions. However, simulations with larger values of S t h is worth, because we 
try to find results which are less dependent on S t h- Moreover, we also try to test whether 
the less-collapsed perturbations are enough to explain Lya absorption. If a model can 
only fit to the observation for S t h > 4.5, it probably implies that the less-collapsed clouds 
are not enough to explain Lya absorption, and a part of the clouds should be located in 
collapsed regions. We should also keep in mind that in the PS formalism the smoothing 
scale can be chosen arbitrarily, but for the baryonic matter considered in our model the 
smoothing scale must be larger than the Jeans length since below the scale thermal pres- 
sure overcomes gravity. Therefore, if we smooth just at the Jeans scale, a threshold value 
might be slightly higher than that in PS formalism. 

Figs. 2c - 7c shows that the model of CHDM always predicts too few absorbers at 
high redshifts for all parameter ranges considered. Even for J-21 = 0.3 and S t h = 7.5 the 
number of the Lya absorbers given by CHDM are still significantly less than the real data. 
In fact, the curves of dN/dz in the CHDM model are almost independent of S t h- This 
shows that the deficit of absorption in the CHDM model is due to the late formation of 
structures. Because the range of S t h > 1.5 is located in the high density contrast tail of the 
perturbation, the thresholds S t h remove only a very small fraction from the original Sjj. 



12 



v(z) is the peculiar velocity in the direction ol the line-of-sight. To synthesize an absorp- 
tion spectrum over a large redshift range, say 1.5 to 4.1, we divide this range into 13 bins, 
each has Az = 0.2 and is centered at z n = n X 0.2 + 1.3 with n = 1 to 13. For a given z n , 
we can use the approximation 5b(x(z), z) ~ 5b(x(z), z n ), because the exponential factor 
in eq.(13) restricts that the integrand has significant contribution to ((Z) only at z ~ Z. 
For each z n , we generate 20 or, sometimes, 100 samples. We take the pixel size in z space 
to be about the same as observations, say, 0.0593A. 

As examples, three simulated absorption spectra at redshift 3 in SCDM, LCDM and 
CHDM are shown, respectively, in Figs. la, b and c. The observational instrumental point- 
spread-function 35 km s -1 is assumed, and the Poisson and background noises are also 
added. Fig. la is for SCDM with J_ 21 = 0.3; Fig. lb for LCDM, J_ 21 = 1.0, and Fig. 
lc for CHDM, J -2i = 0.3. In all three spectra, S t h is taken to be 1.69. The spectra of 
LCDM and SCDM, indeed, appear the features of Lya forest. We can then identify the 
absorption lines from the spectrum by the same way as observation. 

3. Result 

3.1 Number density of Lya clouds and its evolution 

The ^-dependence of the number density of Lya absorption lines with rest equivalent 
width W greater than a threshold Wth is generally described by a power law 



where (dN/dz) is the number density extrapolated to zero redshift, and 7 the index of 
evolution. If the hydrogen distribution is comoving in a flat universe, the number density 
should be (dN/dz) oc (1 + z) 2 /[tt(l + zf + A] 1 / 2 . The deviation of dN/dz from the 
comoving curve implies an evolution of the population of Lya clouds. The index 7 is 
found to depend on Wth- the larger Wth, the higher 7 (Lu et al. 1991; Bechtold 1994). 
Instead of the phenomenological eq.(14), we directly use the observational data of dN/dz 




(14) 
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s -1 Hz -1 sr -1 (e.g. Lu, Wolfe & Turnshek 1991; Bechtold 1994). Recently, J_ 2 i is found 
to be ~ 0.3 at redshift z ~ 4.2 (Williger et al. 1994), which is less than the estimates at 
z ~ 2.5 by a factor of 3 - 10. Therefore, the UV background probably declined from z ~ 2 
to z ~ 4. In our simulation, we will consider several possible models for the background 
flux. 

2.3 Lya forest absorption spectrum 

The profile of the Lya absorption in the spectrum of a QSO at redshift z\ is calculated 
by exp {— t(z/)}, where the optical depth r is given by (Weinberg 1972; BBC) 

t{v) = f° n H i{t)(T a {-)dt , (9) 

where to denotes the present time, t\ the time corresponding to redshift z\ } v the observed 
frequency, which is related to a redshift of Z = (y a jv) — 1, and v a is the Lya frequency 
in rest. The absorption cross section o a is given by 

*« = j-^=V(a,^p), (10) 

0a/7T bv a 



where b = y2kT/m p /c is velocity dispersion in the unit of c, a = 2ire 2 u a /3m e c 3 b = 
4.8548 x 10" 8 /6, I a = 4.45 x 10" 18 cm -2 and V is the Voigt function. Using Eq.(7), t(v) 
can be rewritten as 

T(z) = ^ncn H (i + z) 2 az) , (11) 

'IT 



where ((Z) given by 

C(Z)= I" dzf[l + 6 b (x(z),z)] 1 A n^T^^) , (12) 

and the function x(z) = (1/H ) Jq [fi(l + z) 3 + A] -1 / 2 ^ is a normalized comoving coordi- 
nate. For absorption lines with width W < 0.7A, one can safely replace the Voigt function 
by an exponential profile. We have then 
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C(Z)= I dzf[l + 6 b (x(z),z)] ^ l e -U*-Z)m+Z)+«W?/* , (13) 

/o ^oV (l + z) 3 + X b 
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dimension power spectra of the density and velocity perturbations are given by 



Jk 



'CO 



Pn(k) = 2tt 



P b (q)qdq , 



(5) 



and 



POO 

P lv (k) = 27rH 2 (t)a 2 (t)k 2 / P b (q)q- 3 dq 



Jk 



(6) 



respectively, where ait) denotes the cosmic scaling factor normalized at the present time, 
i.e. a(to) = 1. Discrete one-dimensional sample at redshift z, 6(k } z) } can then be realized 
from the Gaussian probabilities with variances given by Eq.(5). The corresponding density 
fluctuation 5(x,z) in x space can be obtained by the Fast-Fourier transformation. The 
spatial distribution of baryons at redshift z will finally be found by removing all collapsed 
regions from the original 6b(x } z). In the following, we will use the same notation 6b(x } z) 
for the new distribution. Obviously, 6b(x } z) contains also diffuse gas. 

If the abundance of hydrogen and helium is constant over the perturbed distribution, 
the number density of neutral hydrogen atoms is 



where njj is the number density of hydrogen atoms, nc = pc/m p} Pc = 3Hq/8ttG is the 
present critical density, and 0// = X0°, where X = 0.76 being the cosmic abundance for 
hydrogen. The factor / = nui/nn accounts for the neutral fraction. If the medium is in 
ionization balance, one has 



where a, T, and 4:TtJ c Gh I^h are the rates of radiative recombination, collisional ionization 
and photoionization, respectively. The parameters in eq.(8) are: Gh = 2.54 X I0 8 , and 
J v = J c v c jv, where v c is the Lyman limit frequency and J c the external UV background 
flux at frequency v c (Black 1981). From the proximity effect, the Lyman limit background 
flux is estimated to be J_ 2 i = 0.3 to 3.0, where J_ 2 i is defined by J c = J_ 2 i X 10~ 21 cm -2 




(7) 



a(T) 




a(T) + T + 4:TrJ c G H /n H ' 
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2.2 Distribution of Lya absorbers 

All neutral hydrogen located in the line-of-sight of a given distant QSO will contribute 
to the Lya absorption. Therefore, one can model the QSO Lya forest as the absorption 
of HI in the density field calculated under the approximation developed in last section. 

As mentioned in §1, recent observations of the size and velocity of the Lya clouds at 
high redshifts implied that Lya clouds are probably not virialized, and not located in the 
potential well of collapsed objects. Therefore, to model Lya clouds, we should remove 
the regions corresponding to the collapsed objects from the density field of baryonic 
matter. The collapsed regions at high redshift might be related to well developed objects, 
such as CIV and Mgll metal systems and dampled Lya absorption systems. It was 
shown in the Press-Schechter (PS) formalism (1974; also in the peaks formalism, Bardeens 
et al. 1986) that the collapsed structures can be one-to-one mapped from their linear 
correspondences, i.e. it can approximately be identified as regions, in which the mean 
linear density contrasts are exceeding a threshold S t h- We will use a similar threshold to 
remove the collapsed regions. 

Practically, we first do a convolution of the one-dimensional density field by a spatial 
Gaussian filter on the Jeans scale Xb/\/3, where the factor \/3 comes from the projection 
of 3-D to 1-D. We then label collapsed regions on all the peaks where the maximal density 
contrasts are higher than S t h- As in the PS formalism, the parameter S t h essentially is semi- 
empirical. In our case, the only restriction would be S t h > 1.69. To reduce the uncertainty 
given by 5 t h, we will be mainly interested in the results that do not sensitively depend on 
the choice of S t h- 

Therefore, the spatial distribution of baryons involving the QSO Lya forest can be 
calculated in the following way. To find out a one-dimensional distribution of the matter, 
we project all modes of fluctuations in 3-D k-sp&ce into a given 1-D direction. The one- 
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temperature of the Lya forest has been found, one can assume that the gas is isothermal, 
i.e. 7 = 1. We will take the temperature to be T = 5 X 10 4 K corresponding to the velocity 
dispersion b = 29 km/s (e.g. Chaffee et al. 1986; Carswell et al. 1987; Rauch et al. 1992). 
In this case, the molecular weight fi in eq.(l) is ~ 0.64. 

In the linear evolutionary stage, a good approximate solution for the density contrast 
5b and the peculiar velocity v& of the baryonic matter in k-space has been found to be 
(Fang et al. 1993) 

5 h (kJ) - SDMiKt) 



1 + xlk 2 



2 7 kT(t) 1/2 



(4) 



where 5dm ( = Q C 5 C + &h5h) and Vdm are the density and velocity perturbations of the 
dark matter, respectively. The Jeans length xj, of the matter is given by 
1 

b H 3fi//ra p (l + z) 
Eq.(2) shows that the baryonic matter traces the dark matter distribution on scales larger 
than Xb, and smooths over structures on scales less than Xb- The goodness of the approxi- 
mate eqs.(2) and (3) have been checked by numerical solutions in SCDM and CHDM, and 
the solutions were used to model hot (kT ~ 1 keV) component in the intergalactic medium 
(Fang et al. 1993). It was successfully applied to calculate the cosmic Mach number, the 
Comptonization parameter and the soft X-ray background radiation. Therefore, at the 
linear evolutionary stage the density distribution of the baryonic matter can be described 
as a Gaussian perturbation with the power spectrum Pb(k) = Pdm(^)/(1 ~\~ x 1k 2 ) 2 , where 
-Pdm(^) is the spectrum of the dark matter. In the following simulation Pdm(^) will be 
chosen from 1) Davis et al. (1985) for the SCDM model; 2) Klypin et al. (1993) for 
CHDM (since their formula fails at k > 50 Mpc -1 , we extrapolate it by oc k~ 3 ); and 
3) Efstathiou, Bond & White (1992) for LCDM. All the spectra are normalized by the 
COBE's quadrupole anisotropy Q = 13fiK (Smoot et al. 1992). 
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A = tt h = for SCDM; h = 0.5, A = 0, Sl h = 0.30 and fi c + fi 6 = 0.7 for CHDM; and 
h = 0.75, n c + n b = 0.3, and A = 0.7 for LCDM. 



2.1 Linear density and velocity perturbations of baryonic matter 

The difficulty to study the evolution of the baryonic matter at late universe (z <5) 
is the lack of knowledge on thermodynamical or hydrodynamical processes of the compo- 
nent. Since the first generation of luminous objects formed, the baryonic gas was strongly 
affected by energy and mass exchanges from the objects via processes such as heating, 
cooling, shocking, re-ionization, stellar winds and accretions. The details of these pro- 
cesses have poorly been known because of many uncertain factors related to the formation 
and evolution of QSOs and galaxies, and to their injection of energy and mass into in- 
tergalactic space. However, observations have shown that the statistical distributions of 
the Lya forest, say HI column density Njji and Doppler parameter 6, have only a little 
change in the redshift range 2 < z < 4. It strongly indicates that the Lya clouds do pos- 
sess their own cosmological properties. To study these properties except for the proximity 
effect, it is reasonable to consider only cosmologically averaged contribution of the energy 
exchanging processes, but not necessarily to consider the details of individual luminous 
sources. Among these averaged processes, the most important one is the photo-ionization 
of baryons by a mean background UV radiation, i.e. the Gunn-Peterson effect (Gunn & 
Peterson 1965). Therefore, if our interests are limited to the properties of Lya absorbers 
on large scales, the baryonic component can be treated as a cosmic component weakly 
linked to individual luminous objects. Thus, the baryonic matter could be approximately 
described as a polytropic gas with the equation of state p oc p 1 . We have then 



and the mean temperature of the baryonic matter at time t. Since no evolution of the 





where m p is the mass of proton, v s (t) and T (t) are, respectively, the sound velocity 
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The linear approximation will lead to uncertainties in the confrontation between 
observational and simulated samples, especially when clouds have density contrast larger 
than 1. Nevertheless, linear approximation approach is worth: 1) one can find some 
results which are less dependent on the linear approximation, for instance in the cases 
that the variance of density perturbation is less than 1; 2) considering the features of the 
uncertainties given by linear approximation, one can find some instructive constraints to 
models being considered. 

In §2, we describe the method of the simulation. Considering the shortages in the 
works of BBC and B93, we will improve the simulation in the following ways: 1) using 
the density perturbation spectra normalized by the temperature fluctuation in the cosmic 
background radiation (Smoot et al. 1992; Bennett et al. 1994); 2) using a linear threshold 
to identify the regions of collapse and then to remove these regions from the density 
fields. In §3, we give the simulated results, and conducted a confrontation between real 
and simulated samples in various aspects, including the number density and its redshift 
evolution, the equivalent width distribution, clustering, and the Gunn-Peterson effect. 
Discussion and conclusion are shown in §4 finally. 

2. The Method 

We consider a flat universe. The mass density parameters satisfy + A = 1, where 
= J2i^i} &j (j = c } h } b) are the density parameters for cold, hot dark matter and 
baryons, respectively, and A = A/3H 2 with A the cosmological constant and H the 
Hubble constant. Because we are only interested in density perturbations at redshifts less 
than 5, the mass density of the cosmic background radiation field can be neglected. An 
extreme interval of the present baryonic density given by the Big Bang nucleosynthesis 
is fi;, = (0.006 — 0.03) h~ 2 (Copi, Schramm & Turner 1995). In our simulation, we will 
choose fi;, = 0.1. In order to easily compare our results with other studies of structure 
formation (e.g. Jing & Fang, 1994), the parameters at z = are taken to be h = 0.5, 
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their column density should be equal to or larger than 10 17 cm -2 (Mo, Miralda-Escude & 
Rees 1993). Clouds with large size and low column densities are not completely gravity- 
confined. 

Therefore, the Lya clouds are probably neither virialized nor completely gravity- 
confined, but given by pre-collapsed areas in the baryonic density field. This suggests 
that Lya clouds would still be remaining in the stage of evolution that non-linear process 
is less important. Therefore, linear approximation would be able, at least partially, to 
describe the features of Lya absorption clouds. This might be the reason that the linear 
approximation simulation done by Bi, Borner & Chu (1992 hereafter BBC; Bi 1993 here- 
after B93) successfully realized some statistical properties of Lya absorption spectrum. 
Within a certain range of model parameters they found that the simulated spectrum is 
in a good agreement with the observed distribution of line equivalent width in the range 
0.1 A < W < 0.5 A. The two-point correlation function of the synthetical lines is also 
consistent with the observation. 

In this paper, we will develop the technique of simulating Lya absorption in QSO 
spectra in the linear approximation regime. Our goal is to use Lya forest as a cosmologi- 
cal test to discriminate among some popular dark matter models, including the standard 
cold dark matter (SCDM), the cold plus hot dark matter (CHDM), and the low-density 
flat cold dark matter (LCDM). As we knew, objects at higher redshift can serve a very 
promising test to discriminate among models, which have been found in good agreement 
with the observational data of low redshift galaxy clustering and of the microwave back- 
ground radiation (Jing & Fang 1994; Ma & Bertschinger 1994). One can expect that the 
distribution and evolution of Lya clouds on a redshift range from 2 to 4 would also be 
able to play such important role. However, up to now, it is impossible to realize sam- 
ples of ~ 100 Kpc objects in so large redshift range by either hydrodynamical or N-body 
simulation. This is also a motivation of developing the simulation of Lya forest in linear 
approximation. 
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1. Introduction 



Lya absorption line systems, shortward of Lya emission in QSO spectra, indicate 
intervening absorbers with neutral hydrogen column densities ranging from about 10 13 to 
10 22 cm -2 . The absorbers with low column densities, e.g. from 10 13 to 10 17 cm -2 , are 
usually called Lya forest or Lya clouds. In terms of metallicity, clustering, number density 
and redshift evolution, the Lya clouds demonstrate very different properties from higher 
column density absorption systems. These differences probably come from the difference 
in their origins. It is generally thought that the high column density absorption systems 
are most likely due to halos or disks of foreground galaxies, while the low column density 
absorbers are some kind of less clustered clouds consisting of photoionized intergalactic 
gas {e.g. Blades, Turnshek & Norman 1988; Wolfe 1991; Bajtlik 1992). 

Various models for the Lya forest have been proposed, including pressure-confined 
clouds in hot intergalactic medium (Sargent et al. 1980; Ostriker 1988), self-gravitating 
objects (Bahcall & Salpeter 1965; Black 1981), minihalos in cold dark matter (CDM) 
scenario (Rees 1986; Murakami & Ikeuchi 1993), low mass objects formed from the short- 
wave tail in the power spectrum of the density perturbations in the CDM model (Bond, 
Szalay & Silk (1988). These models generally assumed that the Lya forest absorber 
is very well developed object, such as static self-gravitating cloud, virialized baryonic 
cloud in potential well of dark matter, and stably- confined gas in pressure equilibrium. 
However, recent measurements have found that the size of the Lya clouds at high redshift 
is unexpected as large as 100 - 200 h -1 kpc (where h is the Hubble constant in unit 
of 100 km s -1 Mpc -1 ), and their velocity dispersion is unexpected as low as ~ 100 km 
s -1 (Bechtold et al. 1994, Dinshaw et al. 1995). These results cannot matched with 
the pictures of pressure equilibrium and virialization. For instance, if the Lya clouds 
with such large size are well gravitational confined, the Press-Schechter theory show that 
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Abstract 



Our goal in this paper is to test some popular dark matter models by Lya forest in 
QSO spectra. Recent observations of the size and velocity of Lya forest clouds have indi- 
cated that the Lya absorption is probably not given by collapsed objects, but pre-collapsed 
regions in the baryonic density field. Therefore, a linear approximation description would 
be able to provide valuable information. We developed a technique to simulate Lya for- 
est as the absorption of such pre-collapsed regions under linear approximation regime. 
The simulated Lya forests in models of the standard cold dark matter (SCDM), the cold 
plus hot dark matter (CHDM), and the low-density flat cold dark matter (LCDM) have 
been confronted with observational features, including 1) the number density of Lya lines 
and its dependencies on redshift and equivalent width; 2) the distribution of equivalent 
widths and its redshift dependence; 3) clustering; and 4) the Gunn-Peterson effect. The 
"standard" CHDM model, i.e. 60% cold and 30% hot dark matters and 10% baryons, is 
found to be difficult to pass the Lya forest test, probably because it produces structures 
too late and favors to form structures on large scales instead of small scale objects like Lya 
clouds. Within a reasonable range of J v , the UV background radiation at high redshift, 
and 5 t h, the threshold of the onset of gravitational collapse of the baryonic matter, the 
LCDM model is consistent with observational data in all above-mentioned four aspects. 
The model of SCDM can also fit with observation, but it requires a smaller J v and a 
higher 5 t h- This suggests that whether a significant part of the Lya forest lines is located 
in the halos of collapsed objects would be crucial to the success of SCDM. 



Key words: Lya forest - cosmology - dark matter 
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